Haldane charge conjecture in one-dimensional multicomponent fermionic cold atoms 
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A Haldane conjecture is revealed for spin-singlet charge modes in 2A'^-component fermionic cold atoms 
loaded into a one-dimensional optical lattice. By means of a low-energy approach and DMRG calculations, we 
show the emergence of gapless and gapped phases depending on the parity of A*' for attractive interactions at 
half-filling. The analogue of the Haldane phase of the spin-1 Heisenberg chain is stabilized for N = 2 with 
non-local string charge correlation, and pseudo-spin 1/2 edge states. At the heart of this even-odd behavior is 
the existence of a spin-singlet pseudo-spin N/2 operator which governs the low-energy properties of the model 
for attractive interactions and gives rise to the Haldane physics. 
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I. INTRODUCTION 

One of the major advances in the understanding of low- 
dimensional strongly correlated systems has been the so- 
called Haldane conjecture. In 1983, Haldane argued that the 
spin-S* Heisenberg chain displays striking different properties 
depending on the parity of 25. ' While half-integer Heisenberg 
spin chains have a gapless behavior, a finite gap from the sin- 
glet ground state to the first triplet excited states is found when 
2S is even. The Haldane conjecture is now well understood 
and has been confirmed experimentally and numerically. On 
top of the existence of a gap, the spin-1 phase (the Haldane 
phase) has remarkable exotic properties. This phase displays 
non-local string long-range ordering which corresponds to the 
presence of a hidden Neel antiferromagnetic order. ^ One of 
the most remarkable properties of the Haldane phase is the 
liberation of fractional spin- 1/2 edge states when the chain 
is doped by non-magnetic impurities.^ The possibility of a 
similar hidden order has recently been proposed in a differ- 
ent context, by studying the one-dimensional extended Bose- 
Hubbard model."* 

In this Rapid Communication, we will reveal a Hal- 
dane conjecture for spin-singlet modes in a 2Af -component 
fermionic chain at half-filling and for attractive interactions, 
with the emergence of gapless and gapped phases depending 
on the parity of A^. The analog of the Haldane phase is sta- 
bilized for even N with all its well-known properties, while a 
gapless behavior occurs when N is odd. The Haldane physics 
with the alternating gapped or gapless behavior thus translates 
here directly into an insulating or metallic behavior depending 
on the parity of N. To exhibit this even-odd scenario, we will 
consider cold fermionic atoms with half-integer hyperfine spin 
F = N — 1/2 at half-filling (N atoms per site) loaded into a 
one-dimensional optical lattice. Due to Pauli's principle, low- 
energy s-wave scattering processes of spin-F fermionic atoms 
are allowed in the even total spin J — 0,2, ... , 2N — 2 chan- 
nels, so that the effective Hamiltonian with contact interac- 
tions reads as follows^: 
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where ^ is the fermion creation operator corresponding to 
the 2N hyperfine states (a = 1, . • . , 2N) at the site of the 
optical lattice. The pairing operators in Eq. ([T]) are defined 
through the Clebsch-Gordan coefficients for spin-F fermions: 
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case, there are N couplings constants Uj in model ([T]i which 
are related to the N possible two-body scattering lengths of 
the problem. In the following, we will consider a simplified 
version of model ([T]| for N > 2 to reveal explicitly the Hal- 
dane charge conjecture. By fine-tuning the different scattering 
lengths in channel J > 2, we will investigate model ([T]) with 

U2 = ■■■ = U2N-2- 
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and m = J2a'^a,i = 
iCa.i is the density at site i. In Eq. the singlet 
BCS pairing operator for spin-F fermions is y/2N F^q ^ ~ 



= 2U2, V = Uo - U2, 
Za.i is the density at site i. 



^a.i^2N^ 



j, where the 

matrix J' is a 2N x 2N antisymmetric matrix with J'^ = — /. 
When V ^ (Uo ^ U2), model ^ is nothing but the Hub- 
bard model for 2A^-component fermions with a U(2A^)= U(l) 
X SU{2N) invariance. This symmetry is broken down to 
U(l) X Sp(2A^) when V ^ 0.''-^ In the special = 2 case, 
i.e. F = 3/2, there is no fine-tuning and models ([T]l and 
have an exact U(l) x SO(5) symmetry [Sp(4) - SO(5)],^ 
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The zero-temperature phase diagram of model (O away from 
half-filHng has been recently investigated by means of a low- 
energy approach^ '" and large scale numerical calculationsii 
for F = 3/2. In this respect, the physics of F > 1/2 fermions 
is richer than in the standard spin- 1/2 Hubbard chain. '-i^ In 
particular, for U < V < and at sufficiently low density, 
the leading superconducting instability is of a molecular type 
with charge 2iVe.'' " In this Rapid Communication, we will 
show by means of a low-energy approach and density matrix 
renormalization group (DMRG) calculations''* that a Haldane 
conjecture for spin-singlet charge modes emerges in model 
(|2]l at half-filling depending on the parity of N. In the = 1 
case, it is well-known that the half-filled SU(2) Hubbard chain 
displays a critical phase for attractive interaction. The analog 
of the Haldane phase of the spin-1 Heisenberg chain occurs 
for = 2 and attractive interactions. 
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FIG. 1: Phase diagram obtained by the low-energy approach in the 
N = 2 case (see text for definitions); the dotted lines stand for 
(second-order) quantum phase transitions. 



II. STRONG-COUPLING ARGUMENT 

We first give a simple physical explanation of the emer- 
gence of the Haldane conjecture for charge degrees of free- 
dom. It stems from the existence of a pseudo-spin operator 
which carries charge: = ^/nJ2Pq^ ^ and 5f — (jii — 
N)/2. This operator is a Sp(2iV) spin-singlet which is the 
generalization of the Ty-pairing operator introduced by Yang 
for the half-filled spin- 1/2 (i.e. iV = 1) Hubbard model^ 
It is easy to observe that Si satisfies the SU(2) commutation 
relations and generates a higher SU(2) x Sp(2A^) symmetry 
at half-filling along a very special line V = NU. The ex- 
istence of such an extended SU(2) symmetry in the charge 
sector for N = 2 has been first noticed in RefQ. In the 
general N case, one simple way to observe the emergence 
of this symmetry for V — NU is to rewrite the Hamilto- 
nian (|2|i in absence of the hopping term = U{N + 1)]: 

n{t = 0) = 2UJ2z{sf-N{N + 2)/Ay On top of the 

Sp(2iV) symmetry, we thus deduce the existence of an ex- 
tended SU(2) symmetry in the charge sector; moreover, for 
a strong attractive U, the pseudo-spin Si is a spin-iV/2 oper- 
ator, which acts on the degenerate low-lying even occupied 
states (5+)'= |0)^ which one sketches here for N — 1,2: 
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The next step of the approach is to derive an effective Hamil- 
tonian in the strong coupling regime \U\ ^ t. At second order 
of perturbation theory and after a simple gauge transformation 
that changes the sign of the transverse exchange, we find a 
spin-A^/2 antiferromagnetic SU(2) Heisenberg chain: HcS = 
JY.^S^ -S^+i with J = 4t^/[N{2N + 1)\U\]. The Haldane 
conjecture for model (|2]l with attractive interactions thus be- 
comes clear within this strong-coupling argument. When we 
deviate from the V = NU line, the SU(2) charge symme- 
try is broken down to U(l) and in the strong-coupling regime 
the lowest correction is a single-ion anisotropy Dj^ii'^f)'^ 



(with D = 2[U — V/N]). The phase diagram of the result- 
ing model for general N is known from the work of Schulz. — 
For even N, on top of the Haldane phase, Neel and large-Z? 
singlet gapful phases appear while gapless (XY) and gapful 
(Ising) phases are stabilized for odd N in the vicinity of the 
SU(2) line. 

We now turn to low-energy and numerical approaches to 
investigate the strong-weak coupling cross-over and the deter- 
mination of the physical properties of the phases in the vicin- 
ity of the y = 7VC/ line. 



III. LOW-ENERGY APPROACH 



We study here the low-energy approach in the simplest 
F = 3/2 case with the emergence of the striking properties of 
a Haldane insulating (HI) phase. The general N case is highly 
technical and will be presented elsewhere. The low-energy 
procedure for F — 3/2 cold fermions has already been pre- 
sented away from half-filling j^'-S i 'i- In the half-filled case, in 
sharp contrast to the F = 1/2 case, there is no spin-charge 
separation for F > 1/2 since an umklapp process couples 
these degrees of freedom.'^ The exact U(l) x SO(5) continu- 
ous symmetry of model (|2]i is hidden in the bosonization de- 
scription. However, it becomes explicit by a refermionization 
procedure as in the two-leg spin ladder,-^ To this end, we in- 
troduce eight right and left moving real (Majorana) fermions 
— 1, ■ ■ • , 8. The two Majorana fermions account 
for the U(l) charge symmetry, the five Majorana fermions 
^i,...,5 generate the SO(5) spin rotational symmetry whereas 
the last one describes an internal discrete Z2 symmetry 
(ci(4),i ^c^4),^, C2(3),i ^ -«C2(3),i) of modcl ©. Within 
this description, the interacting part of the low-energy Hamil- 
tonian for the spin-3/2 model (|2|i at half-filling reads as fol- 
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with gi,2 = -ao ([/ ± y) , = ao {iU + V) , 54 = aoC^, 
and (75 = qq {U + 2V). The zero-temperature phase dia- 
gram of model (O can then be derived by means of a one- 
loop renormalization group (RG) approach. By neglecting the 
velocity anisotropy, we find the one-loop RG equations: 
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FIG. 2: (Color online): (a-b) Order parameters along the V = —2 
line showing the three different phases CDW, HI and RS. String or- 
ders are computed by taking i and j at equal distance from the center 
of the chain, (c) Various charge gaps A^b (see text), extrapolated in 
the thermodynamic limit from data obtained on L = 16, 32 and 64. 
Quantum phase transitions are located in the grey regions. 



where ga = dga/dl, I being the RG time. The resulting 
phase diagram is presented in Fig. [T] As in two-leg elec- 
tronic ladders, there is a special isotropic ray of the RG flow 
where an approximate SO(8) symmetry emerges in the far 
infrared limit. Along the highly symmetric ray ga — 9 
{a — 1, ... ,5), model (O takes the form of the SO(8) Gross- 
Neveu model which is an integrable massive field theory for 
5 > 0. The resulting gapful phase is two-fold degenerate and 
corresponds to a spin-Peierls (SP) ordering, with lattice order 

parameter Osp = J2i,ai-^yi^i,i^a,t+i + H.c]. A second 
massive phase is obtained from this SP phase by performing 

f7,8 , ^7,8 



which is an exact 



a duality transformation, ^£ 

symmetry of Eq. (|3]l if 54.5 — > —.94.5. This duality symmetry 
exchanges a SP phase with a long-ranged charge density-wave 
(CDW) phase with order parameter Ocbw = J2ii~^T^''^i' 
where 6ni = rii — {ui). The quantum phase transition be- 
tween the SP-CDW phases is found to belong to the U(l) 
universality class. There is a second duality symmetry with 
^£ — >■ — which is a symmetry of Eq. (O if 52,5 ^ ^52,5- 
This duality symmetry is non-local in terms of the original 
lattice fermions Ca.i and gives rise to two non-degenerate 
fully gapped phases from SP and CDW phases. As it is seen 
in Fig. [T] a first non-degenerate phase contains the < 
axis. Its physical interpretation is a singlet-pairing phase 
which is the analog of the rung-singlet (RS) phase of the 
two-leg ladder Upon doping, the singlet BCS pairing /oo.i 
has a gapless behavior and becomes the dominant instabil- 
ity." We need to introduce non-local string order parame- 
ters to fully characterize the last non-degenerate phase. In 
this respect, we define two charge string order parameters: 
0j;vcn ^ COS {^J2k<^^nk), C°f = Sn^O^^J'', which are 
even or odd, respectively , under the particle-hole transforma- 
tion drii — ^ —Srii. Within the low-energy approach, we find 
the long-range ordering of odd (even respectively) charge- 
string operator in the second non-degenerate (RS respectively) 
phase. The phase with long-range ordering of 0°'^'^ is a HI 
phase similar to the Haldane phase of the spin-1 chain. In- 
deed, for attractive interactions U,V <Q, on general grounds. 



we expect that the SO(5) spin gap (As) will be the largest 
scale of the problem. At energies lower than A^, one can in- 
tegrate out the SO(5) spin-degrees of freedom and the leading 
part of the effective Hamiltonian (O simplifies as follows: 



T^int — 



a=7 
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which is the well-known Majorana effective field theory of the 
spin-1 XX Z Heisenberg chain with a single-ion anisotropy 
D?^ Along the special line V — 2U, the two masses mc.o are 
equal due to the presence of the extended SU(2) symmetry 
which rotates the three Majorana fermions ^6'^'^. Within the 
spin- 1 terminology, the interpretation of the phases for U < 
of Fig. [T] reads as follows: the CDW phase is the Neel phase, 
the RS phase is the large-D singlet phase and the HI phase 
is the Haldane phase. All the known quantum phase transi- 
tions in the spin-1 problem are consistent with the findings of 
the RG approach of model (O with a U(l) quantum critical- 
ity for the HI-RS transition and an Ising transition between 
the CDW and HI phases. The HI phase of Fig. [T] is charac- 
terized by a string-order 0°'^'^ which reveals its hidden order 
We can also investigate the possible existence of edge states in 
the HI phase by considering a semi-infinite geometry. In that 
case, the low-energy effective Hamiltonian is still given by 



Eq. © with the boundary conditions: (0) = Cr'" (0)- 

The situation at hand is very similar to the low-energy ap- 
proach of the cut two-leg spin ladder.^^ The resulting bound- 
ary model is integrable and three localized Majorana modes 
77 with zero energy inside the gap (midgap states) emerge in 
the HI phase. These three local fermionic modes give rise to 
a local pseudo spin- 1/2 operator thanks to the identity^^: 
y = —i ?7A ?7/2. We thus conclude on the existence of a 
spin-singlet pseudo-spin- 1/2 edge state which is the main sig- 
nature of the HI phase. 



4 




X 



FIG. 3: (Color online) Integrated excess density F{x) showing the 
edge states in the HI phase of the attractive model with N — 2. 



features of the HI phase are the edge states. As discussed 
above, these edge states are in the charge sector so one can 
observe them by adding two particles while staying in the 
5'^ = = sector. In Fig. [3] we plot the integrated "ex- 
cess density" defined as F{x) — J^dy{n{y) — 2) for U > —1 
and F{x) = {—l)^ J'^dy{n{y) — 2) if [/ < — 1 to remove the 
typical CDW oscillations. We find that, in the HI phase, the 
added particles are pinned at the ends of the chains while in 
the RS and CDW phases, this excess lies in the bulk. 

Finally, we discuss the case = 3, i.e. spin-5/2 fermions. 
As shown in Fig.lH the system behaves effectively as a critical 



IV. DMRG CALCULATIONS 

We now carry out numerical calculations, using DMRG, in 
order to validate this conjecture in the N = 2 and = 3 
cases. We keep up to 2000 states and use open bound- 
ary conditions. When N ~ 2, we fix two quantum num- 
bers for the spin part 5^ = J2a i(-)"+^"a,»/2 and = 
+ "2,j — "3,j ~ «4,j)/2 (the ground state lies in the 
— — sector), and the total number of particles 
Nf = 2L. For TV = 2, we set t = 1, = -2 and we 
investigate order parameters showing the existence of the HI 
phase and its extension. In this respect, we define two string 
order correlations: E{\i - j\) = \{exp{inJ2i<k<]^))\ 

^ndD{\^-J\) = |(i|.exp(z^E.<.<, ^)^)\- InFig.El 
we plot these string correlations, the charge order parameter 
{OcYm{L /2)) in the bulk of the chain, and the pseudo-spin 
gaps which are defined by: 

Aah = E^iNf = 2L + 2b) - E^iNf = 2L + 2a) , 

with Eo{Nf) as the ground-state energy with Nf particles and 
= = 0. In the HI phase, because of the existence of 
edge states (see below), the excited state with Nf + 2 fermions 
falls onto the ground-state (i.e. Aqi = 0), so that the cor- 
rect value for the gap in the bulk is given by A12 ~ A02, 
similarly to what has been done for spin-1 chains. All these 
quantities lead to the conclusion of the existence of two gapful 
phases on top of the CDW phase. In particular, the data con- 
firm the existence of the HI gapped phase with (Ocdw) = 0, 
D{oo) 7^ while E{oo) scales to zero. On the contrary, 
D{oo) = in the RS phase while E{oo) remains finite. In 
the CDW phase, both string orders are finite, which can be 
easily understood from the ground-state structure with alter- 
nating empty and fully occupied sites. One of the striking 




FIG. 4: (Color online) Correlation functions |A''(a;)| and P{x) for 
= 3, corresponding to the correlations of the and S'^ pseudo- 
spin operators, as a function of distance x alV — —6 for L — 72. 
For V = 3U (a), SU(2) symmetry is manifest, while the system is 
XY-like ifU > V/3 (b), or Ising-like ifU < V/3 (c). 

spin-3/2 SU(2) chain on the line V = 3U, with equal trans- 
verse and longitudinal pseudo-spin coiTelations given by the 
singlet-pairing P{x) = (^'oo,L/2+a;^oox/2> and the charge 
coiTelations iV(a;) = ((5ni/2+a;'5'^L/2)i respectively. In par- 
ticular, we recover the same quantum critical behavior as a 
spin- 1/2 chain as predicted.' Moving away from this line, 
we find in Fig. |4]the emergence of a Luttinger liquid phase 
with critical exponents close to the one of the XY model for 
U > V/3, and a gapped Ising phase with exponentially de- 
caying correlations when U < V/3, in full agreement with 
the strong-coupling approach. 
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